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Abstract 

We extend the consideration of the spin- 7; transverse XY chain with correlated 
Lorentzian disorder (Phys. Rev. B 55, 14298 (1997)) for the case of additional Dzya- 
loshinskii-Moriya interspin interaction. It is shown how the averaged density of states 
can be calculated exactly. Results are presented for the density of states and the 
. transverse magnetization. 
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Much work has been done since the famous paper by Lieb, Schultz and Mattis |IJ to 
derive exact results for thermodynamics and spin correlations of one- dimensional spin-| 
XY models. Much less exact results were obtained for random versions of spin-| XY 
chains. One can mention here a group of papers dealing with random spin-| XY models 
using well-known Dyson's and Lloyd's models of disorder [§, £§. Recently the interest in 
random spin-| XY chains has been noticebly increased since they provide a laboratory for 
investigation of generic features of quantum phase transitions in disordered systems. As an 
example we refer to papers on renormalization group and numerical studies on random 
spin-i transverse Ising chain. 

In the present paper we continue the study started in Ref. J| that concerns the spin-| 
isotropic XY chain with random Lorentzian exchange coupling J n and a transverse field Q n 
that depends linearly on the surrounding exchange couplings J n _i and J n . Obviously, due 
to the relation between the transverse field and the random exchange couplings that is a 
model of correlated disorder. The Jordan- Wigner method and the method elaborated 
by John and Schreiber J7| permitted to derive exactly the averaged density of states for 
such a model and as a result to study its thermodynamic properties. Apparently the most 
interesting result of introducing the correlated disorder is the appearance of the nonzero 
averaged transverse magnetization at zero averaged transverse field. Later this effect was 
checked numerically || |§. In the present communication we shall extend the model in- 
troducing additional Dzyaloshinskii-Moriya interspin interaction. Spin-| XY chains with 



Dzyaloshinskii-Moriya interaction were studied in several papers PD| , [TT| , |I2] , |1"3| , [T4|j in which 
it was shown that they exhibit interesting thermodynamic and dynamic properties, which 
may be of interest for the understanding of the properties of some quasi-one-dimensional 
compounds (e.g. CsCuCls). It will be shown below that the Dzyaloshinskii-Moriya inter- 
action may influence in specific manner the thermodynamic properties of a magnetic chain 
conditioned by correlated disorder. 

Hereafter we consider isotropic XY chain in a magnetic field along z axis consisting of 
N spins |. The Hamiltonian is defined by 

N N N 

H = X] ^ nS n + X/ Jn{s n S n+1 + S V n S V n+ i) + ^ D n (s n S y n+1 — S V n S n+1 ), 
n=l n=l n=l 

S n+N = S n- (1) 

Besides the exchange coupling J n between the neighbouring sites n and n + 1 an additional 
Dzyaloshinskii-Moriya interaction D n between these sites is introduced, i.e. a more general 
case than in Ref. f| is considered. 

In what follows we consider two models. 
Model (i) - - We assume the Dzyaloshinskii-Moriya interaction to be ordered, i.e. D n = D, 
whereas the exchange couplings J n are independent random Lorentzian variables with the 
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probability distribution 



P( J n) = -tj M2 - r2 - (2) 

The on-site transverse fields are determined by the formula 

&n ~ — -^(Jn-i + Jn — 2Jq) (3) 

where a is real and | a |> 1. Note that after putting D = one obtains the model considered 
inRef. [|. 

Model (ii) - ■ We assume the exchange coupling to be ordered, i.e. J n = J, whereas the D n 
are independent random Lorentzian variables with the probability distribution 

1 r 

p(D n ) = -— (4) 

7T (D n - Do)^ + 

The on-site transverse fields are determined by the formula 

n n - n = | (D n _i + D n - 2 AO (5) 

where a is real and | a |> 1. 

With the help of the Jordan- Wigner transformation the Hamiltonian (1) can be rewritten 
as a Hamiltonian of non-interacting spinless fermions 

H =Y, n n (c + n c n - ~) + ^ Jn + lD " c + Cri+1 _ Jn 2 lDn c n c+ +1 ^) (6) 

with cyclic boundary conditions. We omitted in (6) the boundary term that is not essential 
for the calculation of the thermodynamic properties ||15|| . Let us introduce the retarded 
and advanced temperature double-time Green functions G^ m (t) = ^fi9 (±t)({c n (t), c+ }), 
^nm(^) = j~ J™oo duje~ lujt G^ m (u ± ie) that satisfy the set of equations 



■ G n-l,m( w ± 1£ ) + o G Z+l,m( u ± le 



uj ± ie - fi n )(7^ m (a; ± ie) 

= <W (7) 



2 n— L,m\ / 1 2 

Our task is to evaluate the random-averaged Green functions since they yield the random- 
averaged density of states through the relation 



p(E) = T-l^GUE±ie). (8) 

7T 

Having the independent Lorentzian random variables one may try to perform the random 
averaging of Eq. (7) with the help of contour integrals. However one must know the positions 
of the singularities of the Green functions in the planes of complex random variables. The 
latter information can be derived for the defined models on the basis of the Gershgorin 



criterion 16 
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Consider at first spin model (i) described by Eqs. (1) - (3). Suppose that exchange 
couplings J n (and hence the transverse fields O n ) are complex variables. As it follows from 
(7) the singularities of the matrix G T =|| G^ m (u> ± ie) || are determined by the zeros of the 
determinant of the matrix A ± iB T where A and B T are the Hermitian matrices given by 



f to — Ref2 x 
-iReJi + \D 



iReJi 
u) — Ref2 2 





-±ReJ 2 - \D 



-\ReJ N + \D\ 




V 



-iReJ, 



N 



uJ 



(9) 



Ref2jv / 



and 



B 



/ e^Imfix ^|lmJi 
^ =F7;ImJ7v 











■ ■ ■ e =F n ) 



(10) 



respectively. John and Schreiber noticed that if all eigenvalues of B T are positive then 
det(A ± iB T ) ^ [[?[]. On the other hand for any eigenvalue A of the matrix B T (10) 
the Gershgorin criterion after making use of Eq. (3) guarantees that at least one of the 
inequalities 



e =F - (Im J n _i + ImJ n ) - A 











< - ImJ n _i - 


f ^ ImJ n | , 


\ a > 1, n = 1, . . 


.,N 



11) 



is satisfied. From (11) it immediately follows that the retarded (advanced) Green function 
does not have poles for ImJ n < (ImJ n > 0) if a > 1 and for ImJ n > (ImJ n < 0) if a < — 1. 

F(. . . , Q — iaT, J — ir, . . .) if F(. . . , Q n , J n , . . .) does not 

= F(...,tt + iaT, J + ir, . . .) if 



Noting that F(. . . ,Q n , J n 



have poles in lower half-planes J n and F(. . . , Q n , J n , . . / 
F(. . . ,Q n , J n , . . .) does not have poles in upper half-planes J n one finds the following result 
of averaging the set of equations (7) 



(u - tt ± i | a | T)G^ m (u) 



Jo — \D =)= i sgn(ct) T 



J Q + iD =f i sgn(a) Y 



G r n+l,m( 1 ^) 



(12) 



The obtained equations (12) possess translational symmetry and proceeding further in stan- 
dard manner one obtains 



p{E) 



1 

7T 



VA 2 + B 2 - A 
\ 2(A 2 + B 2 ) : 



A= (E-Q ) 2 + (1- 



| 2 )r 2 - j 2 



D 



B = 2T[\ a\(E-Q ) + sgn(a) J 



(13) 
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Consider now spin model (ii) described by Eqs. (1), (4), (5). Certainly we may repeat 
the whole calculation once more obtaining as a result p{E) for this model. However there 
is a relationship between the models (i) and (ii) that immediately yields thermodynamics of 
the latter model if it is known for the former one. Namely, consider the following rotations 
of spin axes around z axis 



W = exp 



.n(n - 1) 
-l— st 



s n exp 



.7r(n - 1) 2 
i— ^ — st 



(14) 



One immediately finds that the Hamiltonian (1) arises as a result of transformations (14) 
applied to a Hamiltonian of form (1), however, with the exchange couplings D n and the 
Dzyaloshinskii-Moriya interactions — J n . Therefore it becomes evident that the density of 
states (13) after the replacement Jo — > D , D 2 — > J 2 transforms into the density of states 
for the model (ii). Hence it is sufficient in what follows to consider only the spin model (i) 
defined by (1) - (3). 

Let us discuss the obtained density of magnon states (13). It can be straightforwardly 
checked that (13) covers in the particular case D = the result derived in Ref. In the 
limit of diagonal disorder T— >0, | a | V — 7 = const Eq. (13) reproduces the density 
of states for spin-| isotropic XY chain with Dzyaloshinskii-Moriya interaction in a random 
Lorentzian transverse field with the mean value Qo and the width of distribution 7 Jl4| . The 



density of states (13) remains the same after the simultaneous change of signs of Jo and a; 
hereafter we choose J > 0. 

Let us remind how the density of states is influenced by correlated disorder in case of 
D = (for details see Q). For | a |~ 1 the disorder causes a smearing out of mainly 
one edge of the magnon band (which one depends on the sign of a). As a result we have 
J ^ dE~p{E) Jo 00 dEp{E) at fi = that leads to the appearance of a nonzero averaged 
transverse magnetization = — | /f^ dEp(E) tanh 4p at zero averaged transverse field 
Qq. With an increase of | a | the symmetry of the non-random case is recovered, i.e., both 
edges of the magnon band become smeared out in a symmetric way, the numbers of states 
J ^ dEp(E) and J °° dEp(E) at f2 — become equal to each other, and = at fi — 0. 

Figs, la, lb demonstrate the changes in the behaviour of the averaged density of states 
p(E) versus E — Qq for T = 1, a — ±1.01, Jo = 1 for three different strengths of the 
Dzyaloshinskii-Moriya interaction D = 0, D = 1, D = 2. It can be seen that an additional 
Dzyaloshinskii-Moriya interspin interaction 1) increases the width of the smoothed magnon 
band; 2) leads to the recovering of the symmetry with respect to the change E — Qq —>■ 
— (E — Q ). Thus the increase of the Dzyaloshinskii-Moriya interaction leads to the decrease 
of the nonzero value of WT^ at Sl = (Figs, lc, Id). 

In Fig. 2 we depicted the influence of an increase of the averaged exchange coupling J 
at fixed D = 0. Similarly to the previous case one observes an increasing of the band width, 
however, in contrast to the previous case the density of states remains not symmetric with 
respect to the change E — Q — > —(E — fl ) (Figs- 2a, 2b) and as a result the model exhibits 
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a noticeable nonzero value of at f2 = (Figs. 2c, 2d). The difference in the behaviour 
of the density of states with increasing D or J is not surprising since Jq and D enter in a 
different way into (13). 

To summarize, we have studied the spin-| transverse isotropic XY chain in the presence 
of correlated Lorentzian disorder. Going beyond the results given in Ref. |f] we include 
in the model the Dzyaloshinskii-Moriya interaction. The asumption of correlated disorder 
allows the exact calculation of the averaged density of states p(E). The exact formula 
(13) for p(E) is the main result of the paper. Based on this formula one can calculate in 
a simple way exactly the thermodynamic properties like entropy, specific heat, transverse 
magnetization and static transverse linear susceptibility (see for details H). In that sense the 
presented random quantum spin model may serve as a reference model to study the interplay 
of disorder and quantum effects. In particular, it may be used to test approximations and/or 
calculations for finite systems. As an example we present results for the density of states 
and the transverse magnetization. In particular, we find that the Dzyaloshinskii-Moriya 
interaction may lead to a decrease of the nonzero averaged transverse magnetization at zero 



averaged transverse field that appears due to correlated disorder. It is known |K], [II], |T|, 
13" that in the non-random case the Dzyaloshinskii-Moriya interaction leads to spectacular 
changes in the spin correlations and their dynamics. However, the rigorous consideration 
of correlated disorder in this paper is restricted to thermodynamic quantities based on the 
density of states. The effect of the Dzyaloshinskii-Moriya interaction on the spin correlations 
and their dynamics in the presence of correlated disorder may be studied numerically ]17j] . 
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was presented at the XXth IUPAP International Conference on Statistical Physics (Paris, 
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Conference. He also thanks the University of Magdeburg for hospitality in the autumn of 
1998 when the paper was finished. 



6 



References 

[1] E. Lieb, T.Schultz, and D. Mattis, Ann. Phys. (N.Y.) 16, 407 (1961). 

[2] E. R. Smith, J. Phys. C 3, 1419 (1970). 

[3] H. Nishimori, Phys. Lett. A 100, 239 (1984). 

[4] O. Derzhko and J. Richter, Phys. Rev. B 55, 14298 (1997), and references therein. 

[5] D. S. Fisher, Phys. Rev. Lett. 69, 534 (1992); 
D. S. Fisher, Phys. Rev. B 51, 6411 (1995). 

[6] A. P. Young and H. Rieger, Phys. Rev. B 53, 8486 (1996). 

[7] W. John and J. Schreiber, Phys. Status Solidi B 66, 193 (1974); 
J. Richter, Phys. Status Solidi B 87, K89 (1978); 

K. Handrich and S. Kobe, Amorphe Ferro- und Ferrimagnetika (Academie-Verlag, 
Berlin, 1980) (in German). 

[8] L. L. Gongalves and A. P. Vieira, J. Magn. Magn. Mater. 177-181, 79 (1998). 

[9] O. Derzhko and T. Krokhmalskii, J. Phys. Stud. (L'viv) 2, 263 (1998). 

[10] V. M. Kontorovich and V. M. Tsukernik, Zh. Eksp. Teor. Fiz. 52, 1446 (1967) (in 
Russian) . 

[11] Th. J. Siskens, H. W. Capel, and K. J. F. Gaemers, Physica A 79, 259 (1975); 
Th. J. Siskens and H. W. Capel, Physica A 79, 296 (1975). 

[12] O. Derzhko and A. Moina, Ferroelectrics 153, 49 (1994); 

O. V. Derzhko and A. Ph. Moina, Cond. Matt. Phys. (L'viv) 3, 3 (1994). 

[13] O. Derzhko, T. Krokhmalskii, and T. Verkholyak, J. Magn. Magn. Mater. 157/158, 
421 (1996). 

[14] O. Derzhko and T. Verkholyak, Phys. Status Solidi B 200, 255 (1997); 
O. V. Derzhko and T. M. Verkholyak, Low Temp. Phys. 23, 733 (1997). 

[15] Th. J. Siskens and P. Mazur, Physica A 71, 560 (1974). 

[16] F. R. Gantmachjer, Tjeorija Matrix (Nauka, Moskwa, 1967) (in Russian); 
R. Bellman, Introduction to Matrix Analysis (McGraw-Hill, New York, 1960). 

[17] O. Derzhko and T. Krokhmalskii, Phys. Status Solidi B 208, 221 (1998). 



7 




Figure 1: The density of states (described by Eq. (13)) (Figs, la, lb) and the transverse 
magnetization —rn~ z versus fi at (3 — 1000 (Figs, lc, Id) at fixed J = 1, Y — 1 and 
a = —1.01 (Figs, la, lc) or a = 1.01 (Figs, lb, Id). The short-dashed curves correspond to 
D = 0, long-dashed curves to D — 1 and the solid curves to D — 2. 
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Figure 2: The density of states (described by Eq. (13)) (Figs. 2a, 2b) and the transverse 
magnetization —rn~ z versus VLq aX (3 = 1000 (Figs. 2c, 2d) at fixed D — 0, V — 1 and 
a = —1.01 (Figs. 2a, 2c) or a = 1.01 (Figs. 2b, 2d). The short-dashed curves correspond to 
J = 1, long-dashed curves to J = 1.5 and the solid curves to J = 2. 
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FIG. I. The density of states (described by Eq. (13)) (Figs, la, lb) and the transverse 
magnetization — versus Q at f3 — 1000 (Figs, lc, Id) at fixed Jo = 1, T = 1 and 
a = —1.01 (Figs, la, lc) or a = 1.01 (Figs, lb, Id). The short-dashed curves correspond to 
D = 0, long-dashed curves to D = 1 and the solid curves to D = 2. 

FIG. 2. The density of states (described by Eq. (13)) (Figs. 2a, 2b) and the transverse 
magnetization — fn^ versus H at /3 = 1000 (Figs. 2c, 2d) at fixed D = 0, T = 1 and 
a = —1.01 (Figs. 2a, 2c) or a = 1.01 (Figs. 2b, 2d). The short-dashed curves correspond to 
Jq = 1, long-dashed curves to J = 1.5 and the solid curves to J = 2. 
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